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ABSTRACT 

The rotation of more than 700 pulsars has been monitored using the 76-m Lovell Tele- 
scope at Jodrell Bank. Here we report on a new search for glitches in the observations, 
revealing 128 new glitches in the rotation of 63 pulsars. Combining these new data 
with those already published we present a database containing 315 glitches in 102 pul- 
sars. The database was used to study the glitch activity among the pulsar population, 
finding that it peaks for pulsars with a characteristic age ^ lOkyr and decreases for 
longer values of Tc, disappearing for objects with Tc > 20Myr. The glitch activity is 
also smaller in the very young pulsars (tc < Ikyr). The cumulative effect of glitches, 
a collection of instantaneous spin up events, acts to reduce the regular long term spin- 
down rate of the star. The percentage of \v\ reversed by glitch activity was found 
to vary between 0.5% and 1.6% for pulsars with spindown rates |j>| between 10"^'* and 
3.2 X 10~^^ Hzs~^, decreasing to less than 0.01% at both higher and lower spindown 
rates. These ratios are interpreted in terms of the amount of superfluid involved in the 
generation of glitches. In this context the activity of the youngest pulsar studied, the 
Crab pulsar may be explained by quake-like activity within the crust. Pulsars with 
low spindown rates seem to exhibit mostly small glitches, matching well the decrease 
of their crustal superfluid. 

Through the analysis of glitch sizes it was found that the particular glitching 
behaviour of PSR J0537— 6910 and the Vela pulsar may be shared by most Vela-like 
pulsars. These objects present most of their glitches with characteristic frequency and 
frequency derivative jumps, occurring at regular intervals of time. Their behaviour is 
different from other glitching pulsars of similar characteristic age. 
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1 INTRODUCTION 

Pulsar timing, the method by which the rotation of pul- 
sars is described, is a high precision discipline. Often, with a 
very simple model it is possible to predict every turn of the 
star over many years, with an accuracy of a few microsec- 
onds, or better. This accuracy makes it possible to detect 
and measure very small perturbations affecting the normal 
rotation of the star, supplying information about processes 
inside and outside the pulsar. Two types of timing irregular- 
ity have been recognised which still remain to be well under- 
stood: timing noise and glitches. Timing noise refers to un- 
expected, thus unmodelled, features in the timing residuals 
relative to a simple slowdown model. It can be described as 
a random wondering of the residuals, sometimes presenting 
a clear quasiperiodic behaviour ( Hobbs et al. [2010 1. Perhaps 
most, if not all, of this has recently been shown to arise in in- 
stabilities in the pulsar magnetosphere which result in steps 
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in the slowdwon rate ( Lyne et al.|2010 |. On the other hand, 
glitches are discrete changes on the pulsar rotation rate, of- 
ten followed by a relaxation. Our knowledge of the physics 



behind glitches has increased fairly slowly (Glampedakis & 
Andersson|2009 1 , probably due to the lack of relevant obser- 



vational input to constrain the physical models. The study 
of glitches is important, as they are one of the very few in- 
stances through which we can study the interior of a neutron 
star and the properties of matter at super nuclear density 
(Baym et al. 19691. They also have proved to play an im- 



portant role in the long-term spin evolution of young pulsars 
( |Lyne et al.|1996] ), and, if the gravitational radiation related 
to glitches was to be detected, relevant information on the 



interior and orientation of the pulsar could be obtained ( van 
[Eysden fc Melatos|2008l ). 

Glitches are rare events of very short duration, seen in 
the data as sudden jumps in rotational frequency (//), nor- 
mally ranging between 10~^ /xHz and ~ 100 /liHz. Following 
a glitch, the pulsar sometimes enters a stage of recovery, in 
which the rotation frequency decays towards the pre-glitch 
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value. These recoveries have been interpreted as a signature 
of the presence of a superfiuid in the interior of the star 
( Baym et ah [1969 1. Ghtches are thought to be the result of 
a rapid transfer of angular momentum between this inner 
superfiuid and the outer crust, to which the neutron star 
magnetosphere is attached and whose radiation we observe. 
The crust is thought to be slowed down by electromagnetic 
torques provided by the magnetosphere, and because the 
frictional forces between these two layers are small, the su- 
perfiuid keeps rotating faster than the crust, with an angular 



velocity lag which is reduced during a glitch (e.g. Anderson 
Itoh|fl975 Alpar et al.|1984 l and which subsequently re- 
covers. No change of pulse profile or radio flux density has 
ever been reported to be associated with a glitch in a normal 
radio pulsar. The coincidence observed between glitches and 
enhancement in X-rays flux in a few magnetars (Woods et al. 



2004 



Israel et al.||2007| |Dib et al.||2007j , and particularly 
in the high magnetic field X-ray pulsar PSR J1846— 0258 



(Kuiper & Hermsen 2009 Livingstone et al. 20101, seems 



to belong purely to very high magnetic field neutron stars. 
Outburst episodes, thought to be caused by magnetic field 



decay (Thompson & Duncan 1995 Thompson & Duncan 



1996), are often detected from these objects, but are not 



always found associated with glitches, as glitches are not 
always related to X-ray fiux variations ( Dib et al.|2008 1 . 

As a result of regular monitoring of a large number of 
pulsars, a few of which glitch relatively often, some trends 
in glitching behaviour have been revealed. If a glitch occurs 
when the angular velocity difference between the two main 
inner components is reduced, then the glitch activity should 
decay with the spindown rate of the pu lsar, |z>| (|McKenna| 



fc Lyne|[T990j ). This has been proven by Lyne et al. (20001, 
who showed that the rate of spin-up due to glitches is pro- 
portional to |i>|, for pulsars with \u\ between ~ 10~^* Hzs~^ 
and ~ 10~^^ Hzs~^. The characteristic age Tc — —vjlv has 
also been used as a parameter to describe the glitch activity 
of pulsars, indicating that activity peaks for Tc ~ 10* yr and 



decreases for older pulsars (McKenna & Lyne 1990 Wang 



et al. 20001. Additionally, it has been observed that very 
young pulsars (tc < 2 x 10^ yr) also have little glitch ac- 



tivity ( Shemar fc Lyne|1996 l, an effect attributable to their 
likely higher internal temperature ( McKenna fc Lyne|1990 1 . 

In this paper we present the results of a new search 
for glitches performed using the Jodrell Bank pulsar tim- 
ing database. Together with 186 glitches that can be found 
in the literature, the 128 new glitches found in this work 
are used to study the glitching behaviour of pulsars. Sec- 
tion 2 describes the data and section 3 explains how these 
were analysed to extract the glitches and their main param- 
eters. In section 4 the large glitch database is presented, and 
then analysed in the next section. Finally, sections 6 and 7 
present the discussion and a summary of the conclusions, 
respectively. 



2 THE JODRELL BANK PULSAR TIMING 
DATABSE 

The Jodrell Bank timing database comprises observations 
of more than 700 hundred pulsars, carried out at Jodrell 
Bank observatory (JBO) since 1978, and it is described by 



been performed with the 76-m Lovell telescope, with some 
complementary observations made using the 30-m Mkll and 
42-ft telescopes. Every pulsar is observed at typical intervals 
of 2 to 10 days in a 64-MIIz band centred on 1404 MHz, using 
an analogue filter-bank. Occasionally, observations were also 
carried out in a band centred at 610 MHz. 



3 CHARACTERISING GLITCHES 

Pulsar timing, the method by which the rotation of pulsars 
is described, is based on the analysis of the times of arrival 
(TOAs) of the pulses from the pulsar at the observatory. 
These TOAs are obtained by matching the observed pulse 
profile with a standard, representative template. TOAs are 
corrected for the motion of the observatory around the Solar 
System barycentre and then compared with a simple slow- 
down model describing the rotation of the star, that predicts 
that the pulse number Af, arriving at time t is given by 



= vn(t - to) + -uoit - tof + -i>Q{t - tof + ■ 



(1) 



Hobbs et al. (20041. In summary, observations have mostly 



where vo, vq and Dq are the rotational frequency and its first 
two derivatives at the epoch Iq. The model will give integer 
numbers if the rotational parameters are exactly correct. 
The fractional part of N multiplied by the period of the 
pulsar is called a timing residual. The frequency and its 
derivatives at to are determined by minimising the timing 
residuals, which in an ideal case are expected to be normally 
distributed around zero. 

The signature of a glitch in a plot of timing residuals 
with epoch is normally characterised by the sudden onset of 
a continuous increase towards negative values, relative to an 
ephemeris based upon preceding data, as can be seen in the 
first panel of Fig. [l] where a relatively small glitch is shown. 
A fit of u, V and v to pre-glitch data describes the rotation of 
the pulsar well, but after the glitch a new set of parameters 
is needed in order to minimise the residuals satisfactorily. 
Panel (b) of the figure shows the residuals obtained if we 
attempt to fit the whole data-set in the plot with only one 
set of parameters. The glitch is visible with a characteristic 
shape, showing a sharp cusp-like feature which suggests that 
a single-epoch event happened at that time. Such a pattern 
is characteristic of glitches, while timing noise is normally 
seen as rounded wave-like features in the residuals, without 
a preferred direction and without single-epoch events. 



3.1 Glitch sizes and recoveries 

Glitches can also be seen in a frequency residuals plot, 
obtained after removing the main slope of the frequency. 
Fig. [ijc) shows the evolution of the frequency residuals 
through a glitch. The glitch is observed to be a sudden 
positive step in frequency, in this case followed by a neg- 
ative change of the slope. The size of the frequency step is 
probably the main way to characterise a glitch. It is nor- 
mally expressed as the fractional quantity /S.v/v, where 
is the difference between the frequency after and before the 
glitch. Detected glitch sizes range between 10~^^ and 10~^. 
Glitches in the Crab pulsar are all smaller than 200 x 10~^ 
while in the Vela pulsar almost all glitches are larger than 
1000 X 10"^ ( |Wong et al.|200l||McCulloch et al.||l990 l. 
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Figure 1. A glitch in tiie data of PSR B0531+21, the Crab pul- 
sar. It occurred around MJD 53067 and had a fractional frequency 
jump of Hp/v = 5.33 ± 0.05 x 10-^; a small glitch, (a) The tim- 
ing residuals relative to a slowdown model with two frequency 
derivatives when fitting data only up to the glitch date, (b) Tim- 
ing residuals after fitting all data in the plot; note that the glitch 
feature is still visible. Both these panels have the same scale, cov- 
ering 500 ms. (c) Frequency residuals, obtained by subtracting the 
main slope given by an average v. (d) The behaviour of i> through 
the glitch. 



Most glitches are followed by an increase in the spin- 
down rate |!>|, which may subsequently recover towards pre- 
glitch values. Panel (d) of Fig. [l] shows the evolution of 
V through a glitch. Large glitches and their recoveries are 
easily visualised in a !>-plot. Recoveries can sometimes be 
modelled using an exponential function with a typical time 
constant of ~ 100 days, plus a longer time-scale term, which 
can either be represented by a second exponential with a 
larger time constant (~ 1,000 days) or by a simple linear 
decay of |i>| ( Shemar fc Lyne||1996 l. The step in frequency 
derivative at a glitch is expressed as the fractional quantity 
Az>/!> and detected values range between 10~* and ~ 1. Our 
ability to measure /S.i> depends strongly on how well sam- 
pled the pulsar rotation is around the glitch. Particularly, if 
a set of exponentials are being fitted, the reliability of AO 
will depend on the interval size of the post-glitch data used 



to fit, and of course on the capacity of the model to describe 
the data (see|Zou et al |p008| , concerning PSR B1737-30, 
and Wong et al. (20011, concerning the Crab pulsar). The 
results presented in this paper do not involve fitting of short- 
term recoveries because their parameters depend so critically 
upon the usually poorly known glitch epoch. 

Frequency and v plots are produced by performing fits 
of V, i> and sometimes v to consecutive overlapping groups 
of TOAs, each group typically covering a 200 days interval. 
To produce a series of v and !> values the time interval is 
generally shifted by 100 days and the fit performed again. 



3.2 Glitch detection 

In our searches, all glitches were detected by visual inspec- 
tion of the phase residuals, relative to a slowdown model 
with a maximum of two frequency derivatives. Any feature 
looking similar to those in the top two panels of Fig. [T] 
was considered as a possible glitch, and explored in detail. 
Medium size and large glitches {Av/v > 50 x 10~^) always 
have a clear signature in the timing residuals; they are ob- 
vious and easy to differentiate from timing noise. 

Small glitches are more difficult to identify. The small- 
est glitch ever detected is the one in the millisecond pulsar 
PSR B1821— 24, with a fractional frequency change of only 
0.0095(1) X 10"^ ([Cognard fc Backer|[2004|. In spite of its 



small size, this glitch was easy to detect, due to the natural 
rotational stability and small errors in the TOAs of a mil- 
lisecond pulsar. In contrast, for pulsars with higher levels of 
timing noise and/or larger errors in the TOAs the amplitude 
of the glitch signature could be smaller than the noise vari- 
ations, making detection more uncertain. The size of the 
smallest detectable glitch depends strongly on the timing 
noise levels of the particular pulsar, the signal-to-noise ratio 
of the TOAs, and also on how often the observations were 
made. If the TOAs are typically separated by a time longer 
than the glitch recovery time-scales, or there is a gap with 
no data, then the detection of a small glitch may not be pos- 
sible, or if it is, it might be difficult to get good estimates of 
its epoch and size. 

3.3 Determination of Glitch parameters 

To estimate the glitch epoch and the size of jumps in Av and 
A!>, ephemerides describing the data immediately before and 
after the glitch were built, by fitting v, v and i>, and setting 
the epoch of both ephemerides near to an approximate glitch 
epoch. The two solutions are then compared, and the epoch 
at which the phases are the same is taken as the glitch epoch. 
For a large glitch (e.g. Avjv > 1000 x 10~®) there could be 
several epochs for which this is true between the TOA's T\ 
and T2 that surround the glitch. In this case the glitch epoch 
was taken to be the average of T\ and T2, and the error was 
estimated as {Ti — ri)/4, which corresponds to l-cr of a 
square distribution between T\ and T2. The steps Av and 
Az> are found by taking the difference of the values from the 
two solutions at the estimated glitch epoch. The post-glitch 
ephemeris was always built intending to describe the pulsar 
rotation immediately after the glitch, ideally covering only 
30 or 50 days. Unfortunately this was not always possible, 
because the sampling on some pulsars is rather poor, and 
longer intervals had to be used in those cases. 
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Table 1. All glitches, previously published and found in this work. For re-analysed glitches the new values are shown, and 
the original references are given in the last column. Errors are given in parentheses on units of the last quoted digit. Visit 
http://www.jb.man.ac.uk/pulsar/glitches.html for an online and up-to-date version of this Table. 
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Table 1 — continued 
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4 RESULTS: A GLITCH DATABASE 

This new search for ghtches has found 128 new ghtches in 
63 pulsars, representing an increase of almost 70% in the 
number of known glitches. Together, previously published 
glitches and those found by this search constitute the largest 
glitch database at the time of writing, comprising a total of 
315 glitches in 102 pulsars, including 14 glitches in 6 mag- 
netars. The database contains all the glitches known to us 
at MJD 55000 (June 2009), hence does not include more re- 
cently published events such as those in Yuan et al. (20101 



Av/v {^o-'] 



and Weltevrede et al. (20101. Epochs, fractional sizes hx/jv 



and /S.v/v, and references for all glitches are listed in Ta- 
bleU] while a histogram of all glitch sizes, highlighting the 
new ones, is given in Fig. |2] 

Most of the new glitches found occurred in the period 
MJD 50500-55000, as earlier data had already been anal- 



ysed by Shemar & Lyne ( 1996 1 and Krawczyk et al. ( 2003 1 



However, we have learned to detect and measure smaller 
glitches (e.g. Janssen & Stappers 20061, making it easier 



to identify them in early data, where the signal to noise 
is typically smaller. As a result, 20 small glitches that had 
not been reported before were found in data prior to MJD 
50500. In addition, two relatively large glitches in the pul- 



^ Visit [http://www.jb.man.ac.uk/p uIsar/glltches.htmll for an on- 
line and up-to-date version of Tablejl] 



0.01 0.1 1 10 100 1000 10000 100000 



^ 30 




log(Ai'/i' [10-S]) 



Figure 2. Histogram of the fractional quantity I^u/v for all 315 
glitches. The new glitches are included on top of the previously 
published ones, using a lighter colour, and the contribution of 
magnetars is plotted using black filling. Errorbars correspond to 
the square root of the number of event per bin. 



sar PSR B1930-I-22, observed one after the other in 1987, 
were also measured and included in the database. The rea- 
son that they were not reported before is that data prior 
and between the two events are poor, only allowing measure- 
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MJD 

Figure 3. Frequency, frequency residuals and v evolution of 
PSR B1930+22 during the last 30 years. Frequency residuals arc 
obtained by removing the main slope observed in the top panel. 
Although there are only a few observations between 1980 and 
1987, and v cannot be properly determined, the first two glitches 
are evident in the frequency plots. An even earlier glitch may 
have occured at the begining of 1976, as suggested by the earli- 
est datapoint taken from [Gullahorn &: Rankin] | |1978| and from 
comments in the same work. 



ments of A:/, but no good measurements of Az>. Moreover, 



Gullahorn & Rankin ( 1978 1 reported on observations of this 



pulsar in the period 1975-1976 and mentioned the loss of co- 
herence in the timing residuals towards February 1976. JBO 
observations of this pulsar, starting in 1978, show relatively 
higher frequency and spindown rate, confirming the possible 
occurence of a glitch. The plot in Fig [3] shows the measure- 
ment by Gullahorn & Rankin ( 1978[ ), after being corrected 



to the currently known position of the pulsar, and all JBO 
data. 

During the search, 83 already published glitches were 
re-analysed, and the values obtained during the process are 
included in Table [l] with references to the original papers 
given in the last column. A comparison between published 
estimates of the glitch frequency fractional sizes and those 
obtained in this work are shown in Fig. [4] For glitches 
smaller than C\v jv = 2 x 10~^ there are a few discrepancies 
in fractional size, but in general there is a good agreement 
with published data. 



< 




(AWv)-r 



100 
[10' 



Figure 4. Plot comparing the fractional frequency changes 
{A.v/u) of the 83 glitches re-measured in the course of this work 
with previously published values. The straight line is y = x, indi- 
cating the position of equal estimations. 



number of small glitches is significant between 10"^" and 
10^*, suggesting that, even though their parameters may be 
poorly constrained, most glitches of similar size have been 
detected. 

The analysis presented in the following sections do not 
depend strongly on the completeness of the sample towards 
small glitches. 



4.2 Visualisation of the new glitches 

Plots showing all new glitches are shown in Figs.|5][6) and [7] 
For glitches with ^vjv > 30x 10~^ double-panelled plots are 
presented in Figs. |5] and [6j showing at the top the frequency 
residuals relative to a simple one-derivative slowdown model 
fitted to pre glitch data, and the evolution of i/ through 
the glitch at the bottom. The first two glitches seen from 
PSR B1930+22 are shown separately in Fig.[3] 

Fig. [T] shows the timing residuals of all glitches with 
Au/iy < 30 X 10~^. For small glitches we only show the tim- 
ing residuals relative to a slowdown model using a maximum 
of two frequency derivatives fitted only to pre glitch data. 
For all plots the horizontal axis has been set such that the 
origin corresponds to the estimated epoch of the glitch. 



4.1 Lower limit on detectability 

A lower limit on the detectability of glitch sizes is not sim- 
ple to determine; it varies from pulsar to pulsar depending 
upon TOA errors and temporal coverage. However, a rough 
estimate could be inferred from the whole collection of de- 
tected glitches. The accuracy of glitch size measurements 
seems poor below Aiy/u ~ 10~*, as can be seen in Fig. |4j 
where the comparison with published measurements is plot- 
ted. Nonetheless, the contribution of this work to the total 



5 ANALYSIS 

Fig. [8] shows the location in the period-period derivative 
space (the P-P diagram) of every pulsar for which a glitch 
has been detected. As can be inferred from the diagram, 
glitches are phenomena which are present in many differ- 
ent populations of neutron stars. Only millisecond pulsars 
appear not to glitch, with the exception of PSR B1821— 24, 
which has a relatively large P and one of the smallest charac- 
teristic ages among all millisecond pulsars (tc — 30 Myr) . In 
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Figure 5. Four of the new largest glitches in the Crab pulsar. Every glitch is shown by plotting the frequency residuals (top) and !> 
(bottom) against time. The time axis is measured in days and day zero corresponds to the glitch epoch, which is indicated in MJD below 
the name of the pulsar in each plot. 
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Figure 8. P-P diagram showing with the symbol "x" all pulsars 
that have been seen to glitch. Lines of constant characteristic age 
and lines of constant magnetic field are shown and labelled. 



general, it can be inferred from the figure that most glitching 
pulsars have a characteristic age Tc less than ~ 10 Myr. 



5.1 Glitching rate and the characteristic age 

It has been noted that glitch activity reduces as pulsars age 
i McKenna fc Lyne|1990l|Lyne et al.|2000[ ). Fig |^ shows the 




T, (kyr) 

Figure 9. Number of glitches per year Ng for individual pulsars 
versus the characteristic age for all pulsars observed for at least 
3 years and with one or more glitch detected. The straight line 
is a linear fit to th 
characteristic age. 



number of glitches per year, Ng, for all pulsars known to 
have glitched and which have been observed for at least 3 
years, versus their characteristic age. The number of glitches 
per year was estimated using the whole observation span, 
when this was known, or the time between the first and last 
detected glitch, when no other information was available. 

The most frequent glitchers, labelled on the plot, appear 
to coincide with the objects with more detected glitches, 
showing that the plot is not completely contaminated by ob- 
jects observed for relatively short time. The observed glitch- 
ing rate is clearly correlated with Tc, decreasing for pulsars 
with larger characteristic ages, confirming the trend seen by 
[McKenna fc Lyne| ( [19901 ) in a much smaller dataset. We note 
that due to selection effects these data are not complete to- 
wards low glitching rates. However, this has no effect over 
those pulsars with large glitching rates. 

The maximum glitching rate, for a given value of Tc, can 
be obtained from the envelope defined by the distribution of 
rates in the plot. To describe the envelope, data were binned 
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Figure 6. New glitches satisfying Av/v > 30 x 10~^. Every glitch is shown by plotting the frequency residuals relative to a linear model 
fitted to pre glitch data (top) and v (bottom) against time. The time axis is measured in days and day zero corresponds to the glitch 
epoch, which is indicated in MJD below the pulsar name in each plot. There are two glitches in the top left plot for PSR J0631-+-1036, 
at MJD 50608.277 and MJD 50729. 
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Figure 6 - continued There are three ghtches in the same plot for PSR J1841-0524, at MJD 53562, MJD 54012.88 and MJD 54503. 
There is no much data for PSR J1845— 0316 around the ghtch at MJD 52128, so no good measurements of ;> are possible for this epoch. 
However, the glitch is easily identified in frequency data. 



every half decade of characteristic ago, and the maximum 
gUtch rate per bin was selected. By fitting a straight line to 
these selected values the slope of the envelope was obtained, 
as shown in the plot. The fit indicates that on average, a 
pulsar with a characteristic age Tc will glitch a maximum 
of (6 ± 2) X times per year (with Tc measured in 

kyr). It should be noted that this simple analysis is directly 
related to the number of glitches observed, and has nothing 
to do with the size of the glitches. Accordingly, failing to 



detect small glitches could affect those results. The effect 
of glitch sizes and their frequency are better studied by the 
integrated ghtch activity, which is introduced in the next 
section. 
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Figure 7. Phase residuals (measured in milliseconds) showing different glitches in several pulsars. The time axis is measured in days 
and day zero corresponds to the glitch epoch, which is indicated in MJD below the pulsar name in each plot. 



14 Espinoza et al. 



72 o 



1705-3423 
54408 



-500 





1 1 . 1 1 1 1 








• 






L 1717-16 


• 


■ 51 169 


Si 

< . . , 1 < . . 



-1000 



1000 




-2000 



-1000 



1737-30 

48158 



-100 



-50 



1737-30 
53473.56 



-100 



1737-30 
54695.19 



-100 



1737-30 
54810.9 



-100 -50 



50 



1737-30 
54928.6 



-100 



-50 



1737-3137 
51559 



-100 



1740-31 
49572 



-1000 



1740+1000 

54747.6 



to 

\ 



1757-24 

53033.25 



-400 -200 200 



-100 



\ 



100 



D CN 



1758-23 
50939 



-100 



100 




1809-173 
54365.8 



\ 



1000-500 500 1000 -1000 -500 



500 




500 




-600-400-200 200 
Days from glitch 



-200 200 400 
Days from glitch 



1832-09 

5031 4 



-200-100 100 



1830-08 
L 47541 .3 



1822-C 
: 52810 



-200 



1830-1135 

52367 



-500 



500 



-1000-500 500 1000 



.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 i_ 






1838-0453 


• 


- 54140 

- 


• - 

■ 


-200-100 


100 




-500 



500 



1847-0130 

53426 




-400 -200 200 
Days from glitch 



-600-400-200 200 
Days from glitch 



Figure 7 - continued 
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Figure 7 - continued 



5.2 Integrated glitch activity 

The cumulative effect of spin-up due to glitches on pulsars, 
measured over several years, can be used to study the glitch 
activity and its relationship with other parameters. To in- 
crease the statistical power, in its estimation we consider a 
large number of pulsars, including those that have not yet 



been seen to glitch. Following Lyne et al. (20001, the glitch 



spin-up rate of a group of pulsars is defined as 



fglitch — 



(2) 



where the double sum runs over every frequency jump Avij 
due to the glitch j on the pulsar i, and the sum in the 



denominator is the accumulated years of observation of all 
the pulsars of the group. To calculate the accumulated total 
time of observation of different groups of pulsars we use a 
sample of 622 pulsars that have been observed for more than 
3 yr at JBO (Table|3| plus all pulsars that have glitched and 
are not observed at JBO (Table [2|. 



In Lyne et al. ( 2000[ ) pulsars were grouped according to 
their value of u; each group covering a semi-decade of fre- 
quency derivative. In Table |4] we have reproduced the same 
result using the new glitch database and a larger sample of 
pulsars. The first column contains to the logarithm of the av- 
erage spindown rate of all pulsars in each semi-decade of fre- 
quency derivative. The second column is the total observing 
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Table 3. J-names and MJD ranges of all pulsars observed at JBO with a span greater than 3 years. 



J-name 


Range (MJD) 


.7-name 


Range {M,7D) 


,7-name 


Range 


(M.7D) 


J-name 


Range (MJD) 


J0014+4746 


45120 


54945 


J0754+3231 


44816 


54947 


.71703- 


3241 


47389 


54942 


J1801- 


2451 


47553 


54945 


JUUJb+DoJU 


53249 


54945 


JU / OO — iozcS 


47133 


54939 


J1705- 


1906 


43587 


54935 


J1758- 


2846 


52827 


54926 


JUUo4 — 


44984 


54942 


JUeSl4+74jy 


46277 


54940 


J1705- 


3423 


49086 


54936 


J1801- 


0367 


46719 


54935 


JUUo4 — Uoo4 


48763 


54939 


JUo^U — ioOU 


46118 


54947 


J1709- 


1640 


41332 


54935 


J1768- 


1931 


52686 


54946 


JUU4U+07iD 


46715 


54945 


JUo^3+Uloy 


44816 


54947 


J1708- 


3426 


49086 


54935 


J1801- 


2304 


46694 


54946 


J0048+3412 


46715 


54945 


J0826+2637 


40264 


54949 


J1711- 


1609 


47161 


54935 


J1801- 


2920 


49409 


54947 


JUUOO+oliY 


46715 


54945 


JUo^o — 3417 


43684 


54948 


J1713+0747 


49987 


54940 


J1758+3030 


49872 


54945 


JUU0O+470D 


46720 


54945 


JUeSo7+UOiU 


44808 


54947 


J1717- 


3425 


47880 


54936 


J1769- 


2922 


49086 


54947 


J0102+6537 


46715 


54945 


J0849+8028 


46715 


54946 


J1720- 


1633 


46718 


54945 


J1769- 


2649 


50758 


54945 


J0108+6608 


47387 


54945 


J0846— 3533 


44819 


54948 


J1720- 


2933 


46117 


54942 


J1803- 


2137 


46270 


54948 


J0108+6905 


46715 


54945 


J0855— 3331 


44818 


54948 


J1720- 


0212 


46237 


54935 


J1803- 


2712 


47903 


54947 


J0108— 1431 


49091 


54942 


J0900 — 3144 


52616 


54944 


J1721- 


1936 


48209 


54945 


J18()l- 


1417 


52571 


54946 


J0117+5914 


46715 


54946 


J09()8— 1739 


44815 


54947 


.71722- 


3207 


47161 


,54945 


J18()l- 


2154 


52586 


54940 


J()134 — 2937 


49034 


54938 


.7092 l + (i254 


46715 


54941 


.71721- 


3532 


47907 


.549.34 


J18()5 + 


0306 


47239 


54942 


J0139 + 5814 


44815 


54946 


J0922 + ()G38 


43587 


54947 


J1728- 


0007 


47580 


54945 


J1804- 


0735 


47972 


54942 


J0137+1654 


52754 


54945 


J0927+23 


50059 


54381 


J1726- 


3630 


50682 


54946 


J1802- 


1745 


51610 


52677 


J0141+6009 


45109 


54946 


J0943+1631 


43955 


54941 


J1730- 


3360 


47880 


54946 


J1807- 


0847 


43871 


54942 


JUi47+oyJ2 


46715 


54946 


JUy44 — lo04 


43694 


54942 


J1733- 


2228 


44817 


54945 


J1806- 


1164 


48766 


54948 


JUiOi — UooO 


46241 


54942 


juy4D+uyoi 


46278 


54947 


J1730- 


2304 


49057 


54945 


J1803- 


1867 


51243 


52676 


Tr\1 CO 1 

JUiOJ — iDo7 


43834 


54942 
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50110 


54941 
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46784 


54944 


J1804- 


2717 
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54948 


J0156+3949 


46785 


54945 


J0947+27 


50069 


54941 


J1735- 


0724 


47157 


54944 


J1807- 


2715 


47156 


54948 


J0157+6212 


44818 


54949 


J0953+0755 


40105 


54947 


J1732- 


1930 


49796 


54945 


J1805- 


1504 


52599 


54946 


J()2[)5 + 6449 


52327 


54947 


J1012 — 2337 


43954 


54947 


J1734- 


2415 


52828 


54931 


J1808- 


2057 


46450 


54945 


JU215 + 6218 


50065 


54949 


.11012 + 5307 


49221 


54941 


J1734- 
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54948 
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54940 


Tn'">10 1 A 
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J1809- 
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52922 
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.11300 + 1240 


48138 


54946 


.71738- 


2955 


50787 


54946 


J1812+0226 


47207 


54936 


J0435 + 27 


50173 


54939 


.11311 — 1228 


46716 


54946 


.71739- 


3023 


50728 


.54946 


J1810- 


2005 


50758 


54939 


J0450— 1248 


44817 


54948 


J1321+8323 


45120 


54946 


J1743- 


0339 


43584 


54935 


J1811+0702 


50070 


54936 


Tf\A AO O y1 n 

JU448 — 274y 


49383 


54948 


J loo J — 3UoJ 


49421 


54946 


J1743- 


1361 


46784 


54929 


J1811- 


2439 


52859 


54931 


JU4oJ — i7oy 


46118 


54948 


J looo — 2ooo 


48777 


54942 


J1743- 


3160 


47880 


54926 


J1813+4013 


47155 


54945 


irxAca A 1 CE^^o 
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J 14oo — ooo\j 


48921 
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50803 


54945 


J0502+4654 


46238 


54946 


J1509+5531 


44937 


54946 


J1740- 


3062 


50761 


54925 


J1812- 


1910 


51250 


54945 


J0459— 0210 


49071 


54948 


J1518+4904 


49798 
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J1745- 
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50803 


54945 
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43955 
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54942 


J1744- 


3130 


50761 
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J1816- 


1729 


46271 


54934 


J0528+22()0 


45010 


54947 


.71543 — 0620 


46237 


54940 


J1744- 


2335 


49086 


,54934 


J1816- 


2650 


44818 


54934 


J0534+2200 


43954 


54947 


.11543 + 0929 


44814 


54945 


J1744- 


1134 


49534 


54933 


J1814- 


0618 
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54948 
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54934 




46241 


54947 


Ti fi'-(i^_L')/i 1 a 
.J 1 D o T.) -|- z 4 1 o 


44819 


54945 


.71751- 


3323 






J1823- 


1115 


47475 


54939 


.70630-2834 


45195 


54948 


.71640+2224 


50023 


54946 


.71751- 


2857 


51973 


54946 


J1822- 


1400 


46302 


54936 


.70631 + 1036 


49994 


54942 


.71645-0317 


40486 


54940 


.71753- 


1914 


52828 


54931 


J1820- 


1529 


51245 


54939 


.706,")3 + 80.Sl 


44817 


5494(i 


.71643-1224 


49057 


54940 


.71756- 


2435 


46564 


54945 


J1823- 


3021A 


48269 


54944 


J0700+641S 


44816 


54946 


J1645+1012 


50181 


54941 


J1754+5201 


46722 


54946 


J1823- 


3021B 


47978 


54934 


J06B9+1414 


43955 


54949 


J1648-3256 


49086 


54945 


J1767- 


2421 


43689 


54945 


J1823- 


3106 


44325 


54936 


J0725-163B 


62882 


54949 


J1651-1709 


47389 


54942 


J1764- 


3610 


52718 


54932 


J1821- 


0266 


52732 


54948 


J0726-2612 


62773 


54942 


J1652-2404 


44818 


54942 


J1765- 


1660 


52860 


54931 


J1823+0650 


44815 


54924 


J0729-1836 


43684 


54949 


J1649+2533 


50028 


54946 


J1765- 


2621 


51243 


54945 


J1824- 


1118 


46612 


54936 


J0737-2202 


50832 


54948 


.71650-1654 


49071 


54942 


.71755- 


25211 


51811 


54945 


J1821 + 1715 


50173 


54928 


J0737-3()39A 


52856 


54940 


.71652 + 2651 


50028 


54946 


.71755- 


2534 


51245 


54945 


J 1824 - 


1945 


44817 


54934 


J0737-3()39B 


52971 


54948 


.71654-2713 


49383 


54942 


.71759- 


2205 


46781 


54945 


J 1824- 


2452 


47054 


54934 


.70738-4042 


51895 


54224 


.71659-1305 


46718 


54942 


.71756- 


2225 


51218 


54945 


J1825+0004 


46785 


54947 


J0742-2822 


44838 


54949 


J1700-3312 


49086 


54946 


J1756- 


2251 


52810 


54944 


J1822+0705 


50028 


54936 


J07B1+1807 


49343 


54947 


J1703-1846 


47167 


54935 


J1800- 


2343 


47207 


51461 


J1822- 


0848 


52749 


54939 


. . . continued on next page 
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Table 3 — continued 



J-name 





^^j^^ 




J-name 


Range (MJD) 




J-nanie 




xtange 


^^j^^ 




J-name 




Jrtange 


^^j^^ 


J1825- 


0935 


45008 


54948 




49872 


54936 


TI QC\(\-i-C\'7A(\ 


53511 


54947 


TI Q6;9_L1 Ai n 
jxyo^-j-x'xxu 


46727 


54938 


J1822- 


1252 


51242 


54945 


J1844— 0310 


51610 


54946 


J1907+0918 


51274 


54939 


J1950+05 


49873 


54937 


J1826- 


1446 


46302 


54939 


J1847— 0402 


44816 


54936 


J1909+0007 


44818 


54936 


J1952+3252 


47029 


54945 


J1823- 


0154 


49086 


54947 


J1845— 0826 


52752 


54931 


J1909+0254 


44816 


54936 


J1954+2923 


44808 


54938 


J1826- 


1131 


46270 


54939 


J1845— 1351 


52749 


54931 


J1910— 0309 


44817 


54938 


J1955+2908 


45471 


54938 


J1826- 


1334 


46302 


54944 


J iO'iO — -LOU 1 


52718 


54931 


TI Q1 nj-fl^^S 


47389 


54936 


jxyoo-j-ouoy 


43960 


54938 


J1824- 


0127 


52608 


54932 


,71848 — 0123 


47388 


54936 


J1907+0345 


51643 


54929 


J1957+2831 


50239 


54938 


J1827- 


0958 


46302 


54939 


.71845 — 0316 


51609 


54942 


.11909+1102 


44816 


54938 


J2002+3217 


46564 


54948 


J1824- 


2233 


52859 


54932 


.71845+0623 


52718 


54932 


J1910+1231 


44820 


54938 


J2002+4050 


46717 


54942 


J1824- 


2328 


52751 


54932 


.71845—0743 


51833 


54940 


J1908+0839 


51251 


54944 


J2002+30 


49872 


54940 


J1824- 


1423 




52676 


J ± O^vJ — ± ± ± ^ 


52599 


54931 




51 139 


54933 


T9^^A-l-'^1 '-t7 


44817 


54946 


J1826- 


1526 


51510 


52619 


J1848— 1952 


44817 


54936 


J1909 + 0912 


51252 


54929 


J2006 — 0807 


43591 


54941 


J1829- 


1751 


44818 


54939 


J1847— 0605 


51508 


52621 


J1909+0616 


51509 


54930 


J2005 — 0020 


49384 


54941 


J1827- 


0750 


52974 


54931 


J1849— 0636 


44817 


54936 


J1909+0751 


53693 


54944 


J2006+3104 


53630 


54949 


J1828- 


2119 


52599 


54931 


J 1846— 0749 


52871 


54932 


J1910+0534 


51508 


52620 


J2010+2845 


53630 


54949 


J1830- 


1059 


46612 


54944 




52135 


54942 


TI Q1 nj-l 9^R 


52871 


54948 


T9ni ^-L^R/i'=; 


46718 


54946 


J1828- 


1101 


51243 


54939 


J1848— 1150 


52608 


54931 


J1912+2104 


44817 


54931 


J2018+2839 


45121 


54938 


J1828- 


1057 


51305 


54939 


J1851+0418 


46785 


54925 


J1913— 0440 


40653 


54936 


J2018+34 


53550 


54946 


J1829+0000 


52752 


54932 


J 1848— 0601 


52609 


54932 


J1911— 1114 


49605 


54948 


J2019+2425 


49258 


54936 


J1832- 


0827 


46271 


54939 


,71848+0604 


52733 


54924 


J1914+1122 


44820 


54938 


J2022+2854 


45108 


54938 


J1832- 


1021 


46271 






46723 


54927 


TI (}i '-1-I-1 inn 


47158 


54938 








J183()- 


0131 


52718 


54948 


.71850 + 1335 


46785 


54936 


.71911 + 1347 


52718 


54931 


J2022 + 5154 


45118 


54946 


J1833- 


0827 


46449 


54944 


,71848 — 1414 


49071 


54936 


.71913 + 0832 


51643 


54939 


J2023 + 5037 


46785 


54946 


J1830- 


1135 


51816 


54945 


,71852 + 0031 


46564 


54936 


J1913 + 0446 


51832 


54939 


J2029+3744 


46716 


54941 


J1834- 


0010 


46785 


54381 


,71849 + 0127 


51302 


54948 


J1913+0904 


53004 


54948 


J2030 + 2228 


43587 


54936 


J1833- 


0338 


46782 


54936 




52608 


54932 


TI Qi "^-Li nnn 

J ±y±OT^ xuuu 




54929 


T9n*^'?-L1 7 


OUUZO 


54944 


J1834- 


0426 


44814 


54936 


J1849— 0317 


51610 


54946 


J1913+1011 


51465 


54935 


J2037+3621 


46720 


54946 


J1831- 


0952 


51302 


54939 


J1851— 0114 


52699 


54932 


J1915+1009 


45279 


54948 


J2038+5319 


47388 


54946 


J1831- 


1223 


51478 


52620 


J1851+0118 


52614 


54948 


J1913+1145 


51138 


54925 


J2040+1657 


52751 


54934 


J1831- 


1329 


51250 


54945 


J1854— 1421 


46194 


54936 


J1915+1606 


46671 


54929 


J2043+2740 


50256 


54947 


J1832+0029 


52887 


54947 


11 a^i nn9Q 


53817 


54948 


Ti Q1 ^-L1 


45117 


54930 


T9ndfi flilOl 

J ^U40 — Urii X 


44819 


54941 


J1836- 


0643 


46564 


54936 


J1854+1050 


46724 


54938 


J1914+0219 


52608 


54948 


J2046+1540 


43587 


54932 


J1832- 


0644 


51251 


54939 


J1852 — 2610 


48777 


54936 


J1916+0951 


43587 


54887 


J2048 — 1616 


44818 


54940 


J1833- 


0559 


51510 


54945 


,71853—0004 


52280 


54944 


,11916 + 1312 


47157 


54937 


J2046 + 5708 


46783 


54946 


J1836- 


0436 


46449 


54936 


,71853+0056 


51250 


54948 


,11915 + 0227 


52599 


54948 


J2051— 0827 


49535 


54941 


J1834- 


0602 


51477 




Ti xi'ifi-i-ni 1 


47577 


54948 


TI (i1 7-1-1 


45108 


54947 






54940 


J1837- 


0653 


46270 


54936 


J1853+0505 


51817 


54948 


J1917+2224 


46716 


54932 


J2055+3630 


46057 


54941 


J1834- 


0731 


51833 


54945 


J1853+0545 


51634 


54948 


J1918+1444 


45108 


54945 


J2108+4441 


44817 


54946 


J1834- 


0742 


51670 


54945 


J1853+1303 


53044 


54939 


J1916+0748 


50173 


54930 


J2108— 3429 


49796 


54436 


J1836- 


1008 


46951 


54936 


J1857+0057 


46727 


54938 


J1919+0021 


46001 


54948 


J2113+2754 


44817 


54949 


J1834+10 


49872 


54936 




53630 


54944 


TI Q91 4-1 <34R 


44816 


54940 


J21 13+4644 


44326 


54937 


J1834- 


1855 


52744 


54931 


J1857+0212 


46447 


54944 


J1920+2650 


46718 


54937 


J2116+1414 


44329 


54934 


J1834- 


0031 


52599 


54932 


jl857_j_0943 


46462 


54938 


J1919+1315 


53630 


54925 


J2124+1407 


47131 


54934 


J1836- 


0924 


51510 


54940 


J1856— 0526 


52718 


54931 


J1921+1419 


43587 


54929 


J2124— 3358 


49142 


54940 


J1836- 


0944 


53525 


54932 


J1859+00 


49873 


54936 


J1921+2153 


45118 


54937 


J2139+00 


49873 


54934 


J1835- 


1020 


51250 


54940 




52407 


54947 


TI Q9n-Li 1 1 n 

J X X X X u 


51508 




T9i/i'^ nv'^n 


48918 


54933 


J1835- 


1106 


49071 


54940 


.71858 + 0215 


51510 


52621 


J1922+2018 


44819 


54937 


J2150 + 5247 


46716 


54949 


J1836- 


1324 


53072 


54931 


.71857+0210 


51510 


54948 


J1922+2110 


47139 


54937 


J2149+6329 


44818 


54949 


J1837- 


0045 


49383 


54936 


,71857+0526 


51811 


54948 


J1921+0812 


53083 


54932 


J2155 — 3118 


44819 


54940 


J1837- 


0559 


51153 


54945 


,71900 — 2600 


45121 


54938 


,11924+2040 


45194 


54931 


J2157+4017 


45943 


54946 


J1837- 


1243 


51250 


54948 


T1 QOI -l-fll '^fi 


46724 


54936 


.71926 + 0431 


44819 


54948 


,12212 + 2933 


46716 


54939 


J1837- 


0604 


51089 


54945 


J1901+0331 


44814 


54938 


J1926+1434 


44820 


54937 


J2219+4754 


45943 


54949 


J1838- 


0453 


51251 


54948 


J1901+0716 


46664 


54938 


J1926+1648 


47030 


54940 


J2222+2923 


49872 


54939 


J1841- 


0426 


46270 


54936 


J1901— 0906 


49086 


54936 


J1927+1852 


47236 


54937 


J2225+6535 


44817 


54948 


J1838+1660 


52732 


54926 


J1900— 0051 


51635 


54948 


J1927+0911 


53153 


54944 


J2227+30 


50024 


54946 


J1839- 


0321 


51250 


54946 




44821 


54938 


TI Q9S-I_1 74R 




54936 


T999Q-l-fi9n'=; 


44818 


54935 


J1842- 


0359 


46270 


54936 


,71902 + 0556 


43587 


54938 


.T1929+00 


49872 


54941 


J2229 + 2643 


49829 


54937 


J1839- 


0905 


51670 


54946 


.71902 + t)615 


44817 


54938 


.71932+1059 


40401 


54937 


J2229 + 6114 


51977 


54946 


J1841+0912 


43871 


54936 


.71903 — 0632 


47391 


54936 


.71932+2020 


47155 


54938 


J2235 + 1506 


49219 


54946 


J1839- 


1238 


52599 


54931 


.71900 + 30 


49872 


54945 


.71932+2220 


44816 


54947 


J2242+6950 


46717 


54937 


J1840+564() 


44817 




TI (ifii 1 


49872 


54938 


.J lyOU^Z'i^l 


48505 


54949 


T99/1 c m n 1 


49072 


54946 


J1840- 


0840 


52988 


54931 


.71900 + 0227 


51510 


52620 


J1931+30 


49873 


54948 


J2257+5909 


44817 


54935 


J1840- 


1207 


52599 


54931 


J1901+0355 


52588 


54930 


J1937+2544 


46786 


54937 


J2305+3100 


44817 


54937 


J1841+0130 


52749 


54932 


J1901+0413 


51811 


54929 


J1939+2134 


45297 


54940 


J2305+4707 


46716 


54935 


J1841- 


0346 


51508 


54939 


J1901+0436 


52699 


54941 


J1939+2449 


47236 


54937 


J2308+6647 


44817 


54935 


J1844- 


0433 


46597 


54936 


J1906-0056 


46786 


54938 


J1941-2602 


43584 


54934 


J2313+4253 


43871 


54946 


J1844- 


0538 


46270 


54936 


J1903-0258 


52718 


54932 


J1938+0650 


50028 


54941 


J2317+2149 


44815 


54947 


J1841- 


0624 


51816 


54939 


J1903+0601 


51643 


54929 


J1938+20 


53550 


54948 


J2317+1439 


49868 


54946 


J1841- 


0157 


51670 


54946 


J1906+0709 


46270 


54929 


J1943-1237 


43587 


54934 


J2321+6024 


45197 


54935 


J1844- 


0244 


46447 


54936 


J1904+0004 


49071 


54938 


J1941+1341 


52599 


54935 


J2322+2057 


49258 


54947 


J1842+0267 


52608 


54932 


J1904-0150 


52752 


54932 


J1944-1760 


43591 


54934 


J2325+6316 


43956 


54935 


J1842+0358 


52732 


54932 


J1906+0641 


46270 


54938 


J1945-0040 


43587 


54941 


J2326+6113 


44817 


54932 


J1845- 


0434 


52909 


54936 


J1904+07 


53550 


54925 


J1946-2913 


43587 


54934 


J2330-2005 


45111 


54946 


J1842- 


0415 


51250 


54946 


J1904+080() 


51832 


54930 


J1946+1805 


45118 


54938 


J2337+6151 


46717 


54946 


J1842- 


0905 


51460 


54946 


J1903+0925 


53050 


54948 


.71946+2535 


53630 


54949 


J2346-0609 


49383 


54946 


J1844+1454 


43584 


54926 


J1904-1224 


49419 


54938 


J1949-2524 


43587 


54934 


J2354+6155 


44817 


54946 


J1843- 


0356 


51250 


54940 


J1906+0616 


51608 


54929 


J1946+2611 


52586 


54939 


J2363+2246 


49965 


54947 


J1843- 


0459 


51508 


52620 


J1906+0400 


51492 


54929 


J1948+3540 


44816 


54948 








J1843- 


1113 


51895 


54946 


J1907+4002 


44817 


54945 


J1947+10 


49873 


54937 
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Table 2. Observation time spans (T) of 22 pulsars not observed 
at JBO that have glitched. There are 6 magnetars among them. 



J-name 


T (yr) 


J-name 


T (yr) 


J0146+6145'' 


6 


J1341- 


-6220 


7.4 


J0537-6910 


7.3 


J1539- 


5626 


3 


J0540-6919 


7.6 


J1614- 


5048 


7.2 


J0633+1746 


27 


J1617- 


-5055 


10 


J0835-4510 


37.5 


J1644- 


-4559 


16 


J1048-5832 


8.3 


J1647- 


4552^= 


0.3 


J1048-5937'' 


12 


J1709- 


4429 


3 


J1105-6107 


5.4 


J1708- 


-4009'' 


8.7 


J1123-6259 


5 


J1841- 


-0456= 


3.7 


J1302-6350 


13 


J1846- 


0258 


9.7 


J1328-4357 


1.3 


J23014 


-5852^ 


5 



" 4U 0142+61 

IE 1048.1-5937 
" CXO J164710.2-455216 

IRXS J170849.0-400910 
" IE 1841-045 
f IE 2259+586 




loglvl (lO-'^Hzs-') 

Figure 10. The mean glitch rate {Ng) versus |r> I (Table [ij. The 

solid line is a linear fit to the data and has the form {Ng) oc 
|^|0.47(4)^ 



time span Ti of the group in years and the next four columns 
are the total number of glitches Ng, the mean glitching rate 
(Ng) = Ng/^Ti, the number of pulsars with a glitch de- 
tected Npg and the total number of pulsars in the group 
Np. The 6"* column shows the cumulative effects of the fre- 
quency jumps, which when divided by ^ Ti gives the glitch 
spin-up rate, shown in the last column. Errors for {Ng) are 
estimated as ^ Ng/ ^ Ti and the errors for the glitch spin- 
up rate as Vgutdi/ \/ Ng. The new results, plotted in Fig, 



11 



are perfectly compatible with those obtained by |Lyne et al.| 



(20001. Thanks to the larger sample, we have been able to 
reduce the size of the errorbars significantly, and the u range 
for which a slope close to 1 was claimed is now better de- 
fined (for li'-isl between 10 and ~ 32,000, where !>-i5 is 
z> in units of 10~^^ Hzs"'^). A straight line with a slope of 
one is included in the plot for comparison. We note that the 
largest value of the spin-up rate are the result of the high 
glitch activity of PSR J0537— 6910 and the Crab pulsar. Due 
to the low number of pulsars with very high spindown rate 
in the sample, these two pulsars dominate completely the 
integrated glitch activity of the corresponding v bins (see 
Table |4]|. 



The low glitch activity of pulsars with small spindown 
rates is caused by small and rare glitches. There are no 
glitches detected for pulsars with li'-isl < 0.5 Hz . Fig. 10 
shows the mean glitching rate {Ng) against |z>| using data 
from Tableji] A hnear fit to that data gives {Ng) oc i!>| 
This suggestst that for the first 3 bins no glitches should 
be expected for the given accumulated observation times 
{^Ti), and only one (1.25) glitch should be expected for 
the bin centred at log{|!>_i5|) = —0.73. Using the same fit, 
the necessary accumulated time to observe one glitch can be 
estimated for each bin, and given an assumed glitch size, sen- 
sible upper limits for the glitch activity in those u bins can be 
estimated. It is important to notice that the bins centred on 
logdiz-isl) values of 0.25 and —0.22, which present the low- 
est I'giitch values, have the longest accumulated observation 
times and also the largest number of observed pulsars. There 
is an increase of more than two orders of magnitude of I'giitch 
from these bins towards larger spindown rates. This change 
may imply important differences between objects with lower 
and higher v values. 

The average glitch size for the pulsars in these two bins 
is 0.0004 /iHz, and adding all jumps together the glitch con- 
tribution is not larger than 0.004 /xHz. By using the esti- 
mated accumulated time necessary to observe at least one 
glitch, and assuming a glitch size of 0.0004 /^Ifz, the glitch 
spin-up rate for the 4 bins with lowest spindown rate were es- 
timated and plotted in Fig. |ll^ using white triangles. In the 
same manner, the glitch spin-up rate was estimated again, 
but now assuming that the glitch had a frequency size of 
0.03 /xHz (which roughly corresponds to the centre of the 
general Ai/ distribution, see § 5.31, and the results are plot- 



ted using black triangles. Finally, the glitch activity obtained 
if a glitch with Ai/ = 1 /iHz was detected in the data of 
pulsars in these 4 bins, using the currently available accu- 
mulated observing times, is plotted with a dotted line. As 
gets shorter towards lower the maximum spin-up 
rate estimates grow alike. 

In Fig. |11| 3 the percentage of v reversed by glitch activ- 
ity is plotted as a function of the slowdown rate. For pulsars 
with low values, for which no glitch has been detected, 
the percentage was calculated using the same glitch spin- 
up rates estimated for the upper plot. However, only glitch 
activity corresponding to one glitch with — 0.0004 /iHz 
fall below 1%. If estimated for a glitch activity produced by 
one glitch with a frequency jump of 0.03 /xHz the amount of 
!> reversed by glitch activity rises up to more than 2%, for 
log{|/>|) = 0.018, and about 25% for the first bin. The fact 
that the percentage of i> reversed by glitches of those pulsars 
exhibiting the largest glitch activity is always less than 2%, 
suggests that low spindown rate objects are likely to present 
smaller or similar ratios. Consequently, the glitch activity 
values estimated for one glitch having Af = 0.0004 /iHz 
(white triangles) might represent realistic upper limits. 

5.2.1 Glitch spin-up rate and the characteristic age 

Following the same procedure, the characteristic age Tc can 
be used to divide the sample of pulsars, instead of />. The 
sample was divided in half decades of Tc and the results are 
shown in Table [5] and in Fig |12[ As expected, due to the 
dependence of the characteristic age on v, the glitch spin-up 
rate decreases towards large values of Tc. The curve, similar 
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Table 4. Mean rates of glitch occurrence, (Ng), and mean glitch spin-up rates of pulsars grouped every half decade of 
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to that depending on u, exhibits a change of slope before 
becoming zero. The maximum spin-up rate is found for pul- 
sars with characteristic ages around 10 kyr and is about 
250 X 10"^^ Hzs-\ Because in this case PSR J0537-6910 
is included together with other pulsars of similar character- 
istic age, the maximum glitch spin-up rate is about 3 times 
smaller than the one obtained when grouping pulsars ac- 
cording to their spindown. The glitch activity of pulsars as 
young as the Crab (tc ~ 1 kyr) appears to be lower than 
that of slightly older objects but still large compared to the 
rest of the sample. 



5.3 Glitch sizes 

5. 3. 1 Frequency jumps 

The distribution of the fractional quantity ^vjv of all de- 
tected glitches, shown in Fig.[2j spans almost 8 decades and 
exhibits two gaussian-like overlapping peaks, suggesting that 
glitch sizes follow a bimodal distribution. However, given the 
difhculties associated with finding a general lower limit for 
glitch detection, the completeness of the lower end of the 
histogram is not clear, and the current left edge of the first 
peak could change substantially if many small glitches were 
found close to the present limit of detectability. 

The same bimodal behaviour is also found in the distri- 
bution of the glitch frequency jumps Ai/, shown in the his- 
tograms in Fig. |13| In this plot, to study the details of the 
distribution, data were binned using two different bin sizes; 
either one half a decade wide, or a quarter of a decade wide. 
It can be seen that the overall distribution is composed of a 
wide component covering almost the whole range of At' val- 
ues plus a very narrow peak, centred around Af ~ 25 /xHz, 
involving the largest jumps. 

The contribution of the 6 magnetars that have been 
reported to glitch is also shown in Fig. |13| with the small 
black bins. As already noticed by Dib et al. ( 2008 1, although 



clustering is only due to their very similar low rotational 
frequencies. 



5.3.2 Frequency derivative jumps 

Although measurements of Az> may not be extremely accu- 
rate, for large glitches they are likely to be fairly well es- 
timated. Most estimates of /S.v/v are positive, as expected 
if glitch recovery is understood as an exponential-like de- 
crease to pre-glitch values. Nonetheless, a small number of 
measurements are negative. Most of these are barely signif- 
icant values or they may be the result of sharp features in 
timing noise ( |Lyne et al.|2010"l ), which have masqueraded as 
glitches. 

In Fig. [l4] the frequency jump A.u is plotted against the 
corresponding jump in frequency derivative A!>. The plot in- 
cludes 288 glitches, using black dots for positive /S.v values. 
It is found that in general the two quantities correlate, in- 
dicating that large frequency jumps are generally accompa- 
nied by large frequency derivative jumps. The same glitches 
forming the narrow component in the Ai/ distribution are 
also clustered in this plot, and also show large i> jumps. 



5.3.3 Identifying pulsars responsible for the large glitches 
clustered in the ^ii~IS.v plot 

To identify the objects responsible for these large glitches, all 
pulsars having undergone at least one glitch with a size con- 
tained in the range 10 to 45 /xHz, and accompanied by a A!> 
jump between —450 x 10~^® and —10 x 10~^^ Hz s~^, were 
selected. These criteria identified 20 pulsars, which have had 
a total of 57 glitches satisfying the above restrictions. The 
selected pulsars are listed in Table [6] where the number of 
glitches satisfying the selection criteria are indicated for each 
pulsar, compared to the total number of glitches observed in 
that particular object. The restrictions in Ajy and Ai> come 
from visual inspection of the 1S.V distribution and the Ai>- 



their fractional sizes are all large, contributing to the height 
of the second peak in the fractional size histogram (Fig. [2|, 
they spread in the Ai/ distribution, indicating that their 



Ai^ plot (Fig. 14 1. Glitches with a positive ii jump were not 
considered in the selections. 

The pulsars selected are highlighted in the P-P diagram 
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Table 5. Glitch spin- up rate of pulsars grouped in different ranges of characteristic age. 
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Table 6. Pulsars selected by their glitch Au and Ai/ jumps, 
and the number of events (n) satisfying the selection criteria (see 
text), contrasted with the total number of glitches (Nx)- 
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J0729-1448 
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J1809-1917 


1/1 


B0833-45 


13/16 


B1823-13 


3/5 


B1046-58 


1/3 


J1846-0258 


1/2 


J1357-6429 


1/1 


B1853-I-01 


1/1 


B1610-50 


1/1 


B1930+22 
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B1706- 44 


1/1 
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1/3 
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3/5 


B2334-I-61 
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in Fig. 



15 



which shows all pulsars with P > 1.65 x 10^^^. 
It can be seen that most selected objects are concentrated 
in an exclusive region of the diagram, close to the Vela pul- 
sar, which is plotted with a black diamond. The two objects 
furthest apart, among the selected pulsars, are the X-ray 
pulsar PSR J1846-0258 (the one with the highest P), and 
the relatively old PSR B0355-I-54 (towards the bottom of 
the plot). The first one has exhibited only one large glitch 
{Au = 19 ± 1 /iHz), which was accompanied by magnetar- 



like burst activity ( Gavriil et al.|p008 Kumar & Safi-Harb 



[2008] ); a glitch of this magnitude was not expected in such 
a young pulsar. PSR B0355-I-54 suffered a large glitch back 
in 1986 ( |Lyne 1987) and only small events have been de- 
tected since then (those reported by [Janssen &: Stappers| 



1 2006[ )). The pulsar with the shortest period is the X-ray 
pulsar PSR J0537— 6910, which exhibits high glitch activity, 
with most of its glitches having large Av and AO jumps. The 
diagram includes 2 lines of constant i> enclosing the pulsars 
for which i'gntch/i' 



0.01 (see Fig 111 



5.3.4 Cumulative distributions of glitch sizes 

Whether an individual pulsar tends to exhibit specific glitch 
sizes can be studied using a cumulative distribution of glitch 
sizes. Fig. [16] shows the cumulative distributions of Af 
for the 6 pulsars with more than 10 detected glitches. In 



these plots we can see the probability of occurrence of any 
glitch size. Vela and PSR J0537— 6910 appear to have an 
almost mono-sized glitching behaviour, while, for example, 
PSR B1737-30 and PSR J0631+1036 undergo glitches of al- 
most any size with similar probability. This is reflected in the 



fits performed by Melatos et al. ( 2008 1, where Vela and PSR 
J0537— 6910 stand out, by producing comparatively small 
exponents for the fitted cumulative distributions. Accord- 
ing to that work, different exponents are expected, as they 
would reflect, among some universal properties, the inter- 
nal temperature of the star. Despite this, they acknowledge 
that Vela and PSR J0537— 6910 exhibit exceptional glitch 
size distributions. 

The Av distributions also show that in general Vela and 
PSR J0537— 6910 undergo glitches with frequency jumps 
about 10 times larger than all glitches in the other 4 pul- 
sars. However, we are unable to know whether we have 
observed these 4 pulsars long enough to securely say that 
they do not undergo larger glitches. Neither Vela nor PSR 
J0537— 6910 are possible to observe at Jodrell Bank, but 
we believe that all jumps larger than 0.1 /xHz and 0.5 fiHz 
have been detected for these two objects, respectively. Vela 
is an extremely bright source at radio wavelengths, result- 
ing in good signal to noise, and the smallest glitch reported 
has Au = 0.1 /iHz ( [Cordes et al.|1988[ ). In t he case of PSR 
J0537— 6910, with the information available (Marshall et al 



2004 Middleditch et al. i2006 



it is possible to infer that all 
frequency jumps larger than 0.5 /iHz have probably been 
detected. These limits ensure that their cumulative distri- 
butions are in fact thin and sharp. 



6 DISCUSSION 

In order to discuss the results in terms of a popular glitch 
model, a brief description of that model is presented in the 
next section. The following sections discuss the results. 



6.1 Physics of a glitch 

There are two main models describing the origin of glitches. 
The first regards glitches as star-quakes, produced by rear- 
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Figure 11. Glitch spin-up rate versus the slowdown rate \v\ (top), 
and the percentage of v reversed by glitch activity versus the 
slowdown rate (bottom). The straight line in the upper plot has 
a slope equal to 1 and is not a fit to the data. For low \i>\ values, 
upper limits for the spin-up rate were estimated assuming that 
one glitch with size 0.03 /iHz (black triangles), or 0.0004 /iHz 
(white triangles) happened during the time necessary to have de- 
tected at least one glitch, according to a glitch rate extrapolated 
from the rest of the population (Fig.jToJ. The dotted line repre- 
sents the spin-up rate if one glitch with I^v = 1 /iHz was in the 
data of the pulsars in each bin. The percentage of v reversed by 
the simulated glitch activity plotted with white triangles in the 
upper plot is plotted in the bottom plot, using the same symbols. 



rangements of an oblate crust, which would be evolving to- 



wards a most spherical shape as the star slows down ( Baym 
et al.|1969 |. Glitches in the Crab pulsar could be explained 
by this model, but the higher glitch activity of Vela goes 
beyond the maximum activity that changes of oblateness 
could produce ( |Alpar et al!]|1996[ |Wang et al.||2000[ ). How- 
ever, glitches produced by rearrangements of the crust have 
not been ruled out completely, and they could still be the 
cause of many of the glitches observed, as they could be the 
trigger for the other model described below. 

The second model considers the inner neutron star su- 
perfluid as a reservoir of angular momentum, which is trans- 
ferred to the crust during rapid events, producing what is 
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Figure 12. The mean glitch spin-up rate of pulsars versus char- 
acteristic age for pulsars grouped in semi-decade ranges of char- 
acteristic age. 
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Figure 13. Histograms of the frequency jumps l\v of all glitches 
known. Two bin sizes are used to reveal the features not visible 
with the larger size. Magnetar glitches are also plotted, using the 
small bin size and the darkest colour. 



observed as a glitch ( [Anderson fc Itoh||1975[ ). A rotating 
superfluid is organised as an array of quantised vortices car- 
rying the angular momentum of the whole superfluid body, 
which is proportional to the area density of vortices. Hence, 
if the superfluid is to be slowed down, vortices would need 
to move apart to decrease their density and account for the 
loss of angular momentum. With nothing to stop them, vor- 
tices would be expelled at the outer edge of the body, and 
the superfluid would slowdown normally. However, in the in- 
ner crust these vortices will find forces acting against their 
outward migration, impeding their normal slowdown. The 
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Figure 14. The frequency jump Ai/ of 288 glitches, versus the 
corresponding frequency derivative jump Ai>. GUtches with a neg- 
ative Ai/ are plotted in white filled circles, and those with a pos- 
itive Ai/ are in black. The dashed rectangle encloses those events 
in the range 10 < Aw < 45 /xHz and satisfying —450 < Au < —10, 
in units of 10~^^ Hzs~^. 



neutron superfluid of a neutron star is surrounded by a 
dense lattice of ions stabilised by neutrons which, due to 
the high densities, are also in the form of a superfluid. For 
the moving vortices it is beneficial, in terms of energy, to 
pass through the nuclei, which implies pinning there un- 
til something brakes the pinning force ( |Alpar et al.||1984[ ). 
In particular, because the pinned superfluid is not slowing 
down, the Magnus force between pinned vortices and the 
rest of the superfluid will increase, as it is proportional to 
the rotation rate difference between them. For any pinned 
vortex, given that pinning forces are finite, the rotational lag 
will reach a critical value, for which the vortex will unpin 
and move outwards. The model developed by |Alpar et al.| 
( 1984 1 includes two dynamically different superfluid compo- 



nents. One of them is continuously pinning and unpinning 
from the lattice, driven by thermal fluctuations (called vor- 
tex creep), or quantum tunnelling (if the temperature is too 
low). Therefore, this component is slowing down continu- 
ously, at the same rate as the crust slows down. The sec- 
ond component, involving only a small portion of the whole 
neutron superfluid body, is the crustal superfluid, composed 
of pinned vortices which will unpin only when the Magnus 
force is able to exceed the pinning forces. Gradients in the 
pinning forces can cause over-densities of pinned vortices, 
defining regions with high density of vortices (trap zones) 
and also free vortex regions. The unpinning of vortices in 
one dense region could cause the liberation of vortices from 
other regions in an avalanche like phenomena, producing a 
sudden angular momentum release, that is compensated by 
a sudden spin-up of the crust, i.e. a glitch. The collective 
unpinning of vortices could also be triggered, among oth- 
ers, by crust rearrangements ( Alpar et al. 11996'; 'Middleditch' 
|et al.|2006"l [Melatos et al.|2008| ) or by an instability related 
to the inertial r modes of a rotating neutron superfluid, as 
recently proposed by Glampedakis & Andersson (20091. In 
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Figure 15. Top of the P-P diagram (P > 1.65 X 10~^^), display- 
ing with clear large diamonds those pulsars selected in the Av-Ai/ 
plane (Fig. |14| |. The Vela pulsar is highlighted using a large black 
diamond. In addition to lines of constant characteristic age and 
magnetic field, two lines of constant i> are drawn, indicating the 
extremes of the v range for which t'gutch/l'^l ~ 0.01. 
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Figure 16. Cumulative distributions of glitch sizes (Ai') for the 
6 pulsars with more than 10 detected glitches. The distributions 
have been normalised, and give the probability that these pulsars 
suffered a glitch less than a specific size. 



[Melatos fc Warszawski| ( |2009[ ) the observed glitch activity of 
some pulsars is satisfactorily reproduced when considering 
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unpinning due to thermal fluctuations and Magnus stress 
only, introduced as a noncritical self-organised process. 

6.2 Glitch activity through the pulsar population 
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As previously described by other works ( 


McKenna & Lyne 


1990 Lyne et al.||2000 Wang et al.||2000 


1, we find that the 



with 



and also 



glitch activity of pulsars correlates well 
with Tc. In the frame of the pinning-unpinning model, the 
first relationship is expected, since for a faster spindown the 
angular velocity lag between the crustal superfluid and the 
rest of the star would increase faster, being able to reach a 
critical value in a shorter time. Hence, provided that vor- 
tices have places to re-pin, higher spindown rates should 
produce higher glitch activities. However, the Crab pulsar, 
which has the largest among the glitching pulsars, does 
not exhibit large glitches and its glitch spin-up rate is con- 
siderably smaller than that of PSR J0537— 6910, which has 
a similar |i>|. (compare the last two rows in Table |4) corre- 
sponding almost purely to PSR J0537— 6910 and the Crab 
pulsar, respectively). 

The different glitch activity of these two pulsars could 
be related to age, as proposed by Alpar et al. ( 1996 1 to 
explain the differences between the Crab and Vela pul- 
sars. Although the characteristic age is not an accurate 
age estimator, the age of the respective supernova rem- 
nants confirm that the Crab pulsar (tc ~ 1 kyr) is younger 
than PSR J0 537-6910 (r^ ^ 5 kyr) and the Vela pu lsar 
(tc ^ 11 kyr) ([Comella et al.|1969| [Marshall et al.|1998||As 



chenbach et al. 19951. In general, very young pulsars like 
the Crab, with Tc < 5 kyr, undergo small or medium sized 
glitches {Av < lO/xHz), and show a glitch activity lower 



than older objects, like the Vela pulsar (Fig, 12 l. Addition- 
ally, their glitch activity seems to infiuence the long term 
spin evolution to a lesser extent; the evolution of v in the 
Crab pulsar, though interrupted by glitches, is almost lin- 
ear, while in the case of Vela glitches interrupt the evolution 
almost completely. Perhaps higher temperatures in younger 
pulsars prevent the glitch mechanism from working as ef- 
ficiently as it does for Vela-like pulsars ( [McKenna fc Lyne] 
1990 Link et al.|1999 1. In terms of vortex unpinning, under 
higher temperatures thermal fluctuations could effectively 
compete against pinning forces, and impede the formation 
of large pinning zones. Moreover, Alpar et al. ( 1996 1 sug- 
gested that through quakes the Crab pulsar is still creating 
vortex depletion zones, which by the time the Crab pulsar 
is as old as Vela, it will behave like Vela, producing larger 
glitches and higher glitch activity. In this model, the current 
glitch activity of the Crab pulsar is driven by star-quakes. 



Similarly, Middleditch et al. ( 2006 1 proposed that in general 



young objects are just creating their first surface cracks, that 
generate vortex depletion and trap zones, which will later be 
able to produce the large glitch activity seen in objects like 
Vela. 

The glitch activity is observed to decay steadily towards 
larger characteristic ages (between 5 kyr and lOMyr), and 
the same decay is also observe d a s a function of |!>|. In the 
' (Fig 



range 10 < |!>-i5| < 10" (Fig. [llja,), the slope of the glitch 
spin-up rate is very close to -f 1, implying a simple propor- 
tionality between glitch activity and spindown rate. In terms 
of vortex unpinning, this may reflect the linear dependence 
of the Magnus force with the velocity lag between the crust 
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Figure 17. The magnitudes of glitch frequency jumps Au versus 
pulsar mean spindown rate for all glitches in Table[T] The straight 
line is a fit to the data and has a slope of 0.71(5). 



and the neutron superfluid, supporting this model as the 
main mechanism producing glitches on these pulsars. In such 
a scenario, lower spindown rates would imply larger times 
between glitches, but not necessarily a decrease of glitch 
sizes, which depend directly on the amount of vortices that 
are unpinned. Consequently, no significant change in glitch 
size behaviour is observed for these pulsars, as inferred from 
the plot in Fig. 17 that shows Au versus |i>|. Moreover, a 



monotonous decrease on the number of glitches per year is 
in fact observed in the plot in Fig. |10| that shows the inte- 
grated glitching rate versus spindown rate. 

While the glitching rate seems to decrease monotoni- 
cally as the spindown rate does, the glitch activity presents 
a point where the general slope changes. This is caused by 
the significantly lower glitch spin-up rate of pulsars with 
l^'-isl < 10, which appear to present an abundance of small 
glitches; these pulsars exhibit a large proportion of jumps 
below 0.001 /iHz, a size that is not generally observed in 



higher spindown rate objects (see Fig. 171. The change in 
slope comes with a significant drop of the percentage of u 
reversed by glitches (i'gutch/li'i), from about 1% to less than 
0.01%, as \i>\ decreases (Fig.jnja). 

The ratio t'giitch/|£'| is closely related to Icsi/I, where 
7csf is the moment of inertia associated with the crustal su- 
perfluid, and / is the moment of inertia involved in the nor- 
mal spindown of the star, corresponding to the crust and 



the main superfluid bulk. Ruderman et al. ( 1998 1 assumes 
both ratios to be proportional, while [Link et al.[ ( |1999| con- 
siders !>giitch/|t'| as the minimum possible value for Icsf/I- 
Hence, the amount of crustal superfluid, or its minimum 
amount, varies between 0.5% and 1.6% for pulsars with 
—32000 < !>-i5 < —10, which corresponds to those forming 



the -|-1 slope part of the spin-up rate curve (Fig. 111. Out- 
side this range, the portion of crustal superfluid decreases 
about 2 orders of magnitude. 
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For low spindown rate pulsars, a smaller amount of 
crustal superfluid would agree well with the existence of 
mostly small glitches; and these could be still produced by 
the same mechanism acting on larger spindown rate pul- 
sars. In principle, such a situation should not produce any 
change on the number of glitches per year, as is in fact ob- 



have not been considered. Vela and PSR J0537— 6910 have 



been characterised as quasiperiodic glitchers (Link et al 



served (see Fig. 10 1. The causes of this significant decrease 



of crustal superfluid may be related to a possible evolution 
of the trap zones, due to changes in temperature or other 
internal physical parameters. Two lines of constant |/>| indi- 
cating the upper and lower limits of the -1-1 slope section of 
the glitch spin-up rate curve are drawn in the P-P diagram 
in Fig. |15[ The line crosses the centre of the main bulk of 
pulsars in the diagram. 

For the Crab pulsar (v-is ~ 3.9 x 10^) a very small 
amount of superfluid involved in its glitch activity may ex- 
plain its relatively low glitch spin-up rate, compared slightly 



higher spindown rate objects (Fig. 111. In opposition to low 



spindown objects, glitches in the Crab pulsar are not partic- 
ularly small, a fact that may support the hypothesis propos- 
ing quake-driven glitch activity on the Crab pulsar and other 
very young neutron stars. 

6.3 Glitch size distributions 

The study of the frequency jump distribution of all de- 
tected glitches showed that most glitches present a frequency 
jump of between 0.001 and ~10 /xHz in size. Additionally, a 
number of events exhibit larger jumps, narrowly distributed 
around ~ 25 /iHz. These glitches also show large Ai> 
jumps, and are found tightly grouped in a Au-Af plot, al- 



most isolated from all other glitches (Fig. 14 1. Most pulsars 



responsible for these large glitches are Vela-like pulsars, in 
the sense that they have similar rotational parameters, im- 
plying similar characteristic ages and magnetic fields. Those 
pulsars satisfying the selection criteria applied to identify 
these objects, but that are not Vela-like pulsars, may be- 
long to the wide component of the overall Av distribution. 

The cumulative distribution of glitch sizes in the left 
plot in Fig. [16] suggested that there may be two styles of 
glitching, leaving Vela and PSR J0537— 6910 as representa- 
tives of a special class, that mainly presents glitches with 
a similar size, having low probability for smaller events. In 
contrast, pulsars hke PSR B1737-30 and PSR J063H-1036 
present broad glitch size distributions. Presumably, many of 
the pulsars selected by their Au and Ai/ values present sharp 
and narrow cumulative glitch size distributions, like Vela and 
PSR J0537-6910. In Table [6] the number of glitches satis- 
fying the selection criteria per pulsar are indicated, along 
with the total number of detected glitches for that partic- 
ular pulsar. Objects like PSR B1757-24, PSR B1800-21, 
PSR B1823-13 and PSR B1930-I-22 seem to have a large 
proportion of large glitches over small ones. On the other 
hand, PSR J0729— 1448 has only had one large glitch, after 
5 very small events. In general, pulsars exhibit very different 
glitching behaviours and it is difficult to associate them ac- 
cording to their glitching properties. However, Vela and PSR 
J0537— 6910 present clear similarities that may be shared 
with other Vela-like pulsars, suggesting they have something 
in common, which is not found in the rest of the population. 

While the above arguments are related to the specific 
values of the frequency jumps, the times between glitches 



'1999; Middlcdi tch et al.[ 2006: Mclatos et al.|2008| ), because 
their glitches occur at semi-regular intervals of time. There 
are indications of the same regularity between large glitches 
in the Vela-like pulsars PSR B1757-24, PSR B1800-21, 
and PSR B1823— 13, with the average time between glitches 



varying from pulsar to pulsar. Mclatos et al. (2008 1 analysed 



the size and waiting times distributions of a number of pul- 
sars and related their glitch activity to avalanche dynamics. 
In this context the existence of quasiperiodic glitchers comes 
natural and it is seen as the result of mean-field forces dom- 
inating local interactions, which in this case would refer to 
the crust slowdown and thermal fluctuations. 

The different glitching properties between Vela-like pul- 
sars and other glitching pulsars with large spindown rate do 
not seem to be related to the amount of crustal superfluid, 
as these objects are well integrated in the general glitch 
spin-up rate trend, defined by most pulsars in the range 
-32000 < !>_i5 < -10 (Fig. [lT|. Vela-like pulsars seem 



to enjoy a very stable situation, in the sense that the un- 
pinning of a very similar amount vortices (i.e. similar glitch 
sizes) is produced at regular intervals of time. Other objects, 
like B1737-30 and J0631-I-1036, even though with a similar 
characteristic age to Vela, appear to be in a more chaotic 
state, where glitches of any size can happen at any time. 
This could be understood as an inhomogeneous distribution 
of pinning zones, with different pinning capacities; or as the 
absence (or small presence) of whichever dynamic dominates 
the glitch activity of Vela-like pulsars. At present, it is not 
clear whether these differences are part of the normal evo- 
lution in the life of pulsars, or these are exclusive properties 
of a different class. 



7 CONCLUSIONS AND SUMMARY 

A new search for glitches in the rotational behaviour of ra- 
dio pulsars, using the Jodrell Bank pulsar timing database, 
found a total of 128 new glitches in 63 pulsars. These glitches 
plus those already published constitute the largest glitch 
database yet assembled, containing 315 glitches in 102 pul- 
sars. 

The glitch database and a sample of 622 pulsars ob- 
served for at least 3 years at JBO, were used to estimate the 
glitch spin-up rate (fc'gutch) as a function of the characteris- 
tic age (tc), and as a function of the first derivative of the 
spin frequency (v). The glitch spin- up rate peaks for pulsars 
with Tc ~ 5 kyr, and as Tc increases, the rate decreases 
linearly (in a logarithmic space) over 4 orders of magnitude 
until Tc ~ 5 Myr. For longer Tc the glitch activity decreases 
with a higher rate and it disappears for Tc > 20 Myr. To- 
wards the other end, for Crab-like pulsars (tc < 5 kyr), the 
glitch activity seems to be lower than for Vela-like pulsars 
(tc ~ 10 kyr), an effect attributable to higher temperatures 
or insufficient vortex trap capacity. 

The amount of v reversed by the cumulative effect of 
glitches varies between 0.5% and 1.6% for pulsars having 
slowdown rates It'-isl between 10 and 32,000, which includes 
all pulsars with a characteristic age between 5 and 100 kyr 
as well as a few other older objects (see the P-P diagram 
in Fig. |15[ where lines of constant p and Tc indicating these 
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limits are drawn). Towards both extremes, for faster and 
slower spindown rates, this ratio decreases quickly, reach- 
ing values around 0.01% (Fig.jTTjj). In the pining-unpinning 
model these percentages may be indicative of the amount 
of superfluid that is involved in the glitch activity. In this 
context, the linear increase of i>giitch with |i>| can be un- 
derstood as a direct consequence of the linear dependence 
of the Magnus force with the velocity lag between the two 
inner components of the star. 

Low spindown rates combined with small portions of 
crustal superfluid could explain the presence of mostly small 
glitches among low |i>| pulsars. On the other hand, the glitch 
activity of high or very young objects like the Crab pul- 
sar, appears well explained by crust rearrangements, quakes 
or crack growing like models. In this way, despite having 
small portions of crustal superfluid, these pulsars could still 
suffer of medium sized glitches and exhibit a relatively low 
glitch spin-up rate, as observed. 

Among the pulsars for which the percentage of the slow- 
down reversed by their glitch activity is around 1%, the 
study of glitch sizes showed that there are at least two dif- 
ferent glitching styles. One exhibits glitches of any size, typ- 
ically exhibiting 0.0001 < < 10 /iHz, at random in- 
tervals of time; PSR B1737-30 and PSR J0631-I-1036 are 
good representatives of this kind of activity. The other one 
is restricted to an almost unique glitch size, between 10 
and ~ 45 /iHz, accompanied with AO jumps in the range 
-450 X 10"^^ < A!> < -10 X 10"^^ Hz s"\ and occurring 
at semi-regular intervals of time. The pulsars responsible for 
this glitch activity are mostly Vela-like pulsars, occupying an 
exclusive place in the P-P diagram. The differences between 
the glitching properties of Vela-like pulsars and some other 
glitching pulsars, like PSR J1737-30 and PSR J0631+1036, 
suggest that either there are different classes of pulsars, or 
there is an evolutionary trend, being a moment in the life of 
pulsars at which the glitch activity turns into Vela-like. In 
this last scenario, the evolution of the spatial distribution 
and relative sizes of vortex trap zones may be of relevance 
to explain the transition between the two different glitching 
styles observed. 
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